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Abstract. We study the iterated quasi-reversibility method to regularize ill- 
posed elliptic and parabolic problems: data completion problems for Poisson’s 
and heat equations. We define an abstract setting to treat both equations at 
once. We demonstrate the convergence of the regularized solution to the exact 
one, and propose a strategy to deal with noise on the data. We present nu¬ 
merical experiments for both problems: a two-dimensional corrosion detection 
problem and the one-dimensional heat equation with lateral data. In both 
cases, the method prove to be efficient even with highly corrupted data. 


1. Introduction. We consider data completion problems for elliptic and parabolic 
operators. We start with elliptic operators: we consider a bounded domain 0 C 
d > 2, with Lipschitz boundary (see [3]). Let v G be the exterior unit 

normal of dfl, and T, Tc C dfl, such that dfl = T U Tc and meas(r), meas(rc) > 0. 
Let (7 : fl I—>■ be a real matrix valued function such that cr G and 

cr = c o,.e. in fl. 

The data completion problem is: 

Problem. For f, go and in L^(C) x L^(r) x L^(r), find u G such that 

{ —V • tjVw = f in n 
u = go onV 

aVu ■ V = gM onT. 

This problem is well-known to be ill-posed (see [1, 2] and the references therein): 
it does not necessarily admit a solution for any data (/,5 _d, <7Ar), and if a solution 
exists, it does not depend continuously on the data. On the other hand, if the 
problem admits a solution Ug, this solution is necessarily unique (see e.g. [1, 4]). 

Such a problem is encountered in many practical applications, among others in 
plasma physic [5, 6], or corrosion detection problems [8, 9, 7, 11, 10]. We will be 
particularly interested in the corrosion detection problem: in this problem, u is the 
electrical potential inside a conductive object 12, a is the conductivity of the object, 
represent a current imposed on T, accessible part of the boundary of 11, and 
go is the corresponding potential measured on T. The aim is to determine if some 
portion of the inaccessible part of the boundary Tc is corroded. 
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Mathematically, it exists a non-negative function /j. define on such that 

a'S/u ■v + ^j,u = 0 on Tc 

and the objective is to reconstruct fi: /i = 0 on the healthy part of Fc, and ^ > 0 
on the corroded part. In section 6.1, we test our method on this problem. 

The data completion problem is known to be severely, even exponentially ill- 
posed [2] . Therefore ones needs to use regularization methods to try to reconstruct 
u. Several methods have been proposed to stabilize the problem: see, e.g., [12, 13, 
14, 15, 16, 17] and the references therein. 

We are also interesting in the data completion problem for the heat equation, 
which is quite similar to the elliptic one, except that this time u solves a parabolic 
equation. Such inverse problem appears naturally in thermal imaging [33] and 
inverse obstacle problems [34, 35]. For T > 0, we define Q := (0,r) x fl. Let / be 
in L'^{Q), go and gN in L^(0, T; L^(F)). The data completion problem is then 

Problem, find u £ := L‘^{0,T-, (1 H^{0,T-such that 

{ dtu — Am = f in Q 

u = go on (0, T) x F 
Vm ■ V = gN on (0, T) x F 

This parabolic data completion problem is also severely ill-posed (see e.g. [18]). 
Note that it is not mandatory to impose an initial condition m(0, .) on 11 to obtain 
the uniqueness of the solution (if such a solution exists). Again, regularization 
methods are needed to obtain a stable reconstruction of u from the data /, go and 
9n- 

The quasi-reversibility method is such a regularization method, introduced in 
the pioneering work of Lattes and Lions [19] to regularize elliptic, parabolic (and 
even hyperbolic) data completion problems. The mean idea of the method is to 
approach the ill-posed data completion problem by a family of well-posed varia¬ 
tional problems of higher order (typically fourth order problems) depending on a 
(small) parameter e. The solution of the regularized problem converges to the so¬ 
lution of the data completion problem, when the parameter e goes to zero. The 
quasi-reversibility method presents interesting features: first of all the variational 
problems appearing in the method are naturally discretized using finite element 
methods^ thus the method can be used in complicated geometries, an interesting 
property when the method is used in an iterative algorithm with changing domain. 
Furthermore, the method is independent of the dimension. Since its introduction, 
the quasi-reversibility method has been successfully used to reconstruct the solu¬ 
tion of elliptic [20, 21, 23, 24, 25] and parabolic [29, 30] ill-posed problems, and as 
a keystone in the resolution of inverse obstacle problems in the exterior approach 
[26, 27, 28]. 

In the present paper, we are interested in a natural extension of the quasi¬ 
reversibility method, the iterated quasi-reversibility method: it consists in solving 
iteratively quasi-reversibility problems, the solution of each one depending on the 
solution of the previous one. We therefore obtain a sequence of quasi-reversibility 
solutions, which converges to the exact solution of the data completion problem if 
exact data are provided, for any choice of the regularization parameter e. This has 
interesting consequences from a numerical point of view: first of all, one can now 
choose a large value for the parameter of regularization e, leading to an improvement 
in the conditioning of the finite-element problems, without lowering the quality of 
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the reconstruction. This is not the case for the standard quasi-reversibility method, 
for which it is mandatory to use small e to obtain a good reconstruction. Further¬ 
more, in presence of noisy data, we present a method to choose when to stop the 
iterations according to the amplitude of noise on the data, based on the Morozov 
discrepancy principle, which ensure both stability and convergence of the method. 
The main drawback of this extension of the quasi-reversibility method, compara¬ 
tively to the standard quasi-reversibility, is that several problems have to be solved 
to obtain a good reconstruction. However, as it is the same variational problem 
that appears in each iteration of the method, one can precompute a factorization of 
the finite-element matrix. Hence, the cost of the method is not significantly higher. 

The paper is organized as follows. In section 2, we introduce an abstract setting 
to treat both data completion problems we are interested in at once. In section 3, 
we present the standard quasi-reversibility regularization in this abstract setting, 
and prove some results we need to study the iterated quasi-reversibility method. In 
section 4, we focus on the iterated quasi-reversibility method, both in the case of 
exact data and noisy data. In section 5, we show that the abstract setting apply 
to both elliptic and parabolic data completion problems. In section 6, numerical 
results are presented, demonstrating the feasibility and efhciency of the method for 
both problems. 

2. An abstract setting for data completion problems. In this section, we 
set up an abstract setting corresponding to both data completion problems we are 
interested in. 

Let X, y be two Hilbert spaces endowed with respective scalar products (., .)x 
and (.,.)y, and corresponding norms denoted ||.||;r and ||.||y. 

Let y G y. Both of our data completion problems can be written in the following 
way: find x € X such that Ax = y, with A : X y a continuous linear operator 
with following properties: 

• A is one-to-one 

• A is not onto 

• 

In this setting, y plays the role of the data, and x the solution of our data completion 
problem. The problem is obviously ill-posed: indeed, as A is not onto, there exist 
y in y for which the problem admits no solution. We define 3^odm := lm(A) the 
set of admissible data, and ynadm = 3^ \ 3^odm the set of non-admissible ones. By 
definition, yadm is dense in y. Actually, this is also true for ynadm 

Proposition 1. The set ynadm is dense in y. 

Proof. This is quite simple: suppose it exists y € yadm and i5 > 0 such that 
\\y - y\\y <5^ye yadm- It exists x e X s.t. Ax = y. 

y — y d _ 

Let y be any element oi y, y A y. We define y = — - -b y. Obviously, 

\\y-y\\y^ 

||y~y||y < d. Therefore, y G yadm, and it exists x & X such that Ax = y. A simple 
computation shows then that 

=y. 

Hence Im{A) = y, contradicting the assumptions on A. Therefore, for any y G 
yadm, for any (5 > 0, there exists y^ G ynadm such that ||y — y^\\y < d, which ends 
the proof, as 3^ = yadm O ynadm- ^ 
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In other word, for any admissible data y exists a non-admissible one y arbitrary 
close to y. In particular, this leads to the high instability of the problem with 
respect to noise: 


Proposition 2. For any y € y, the exists a sequence Xn € At such that 

n—>-oo 


la^nllA +00 and Axn 


y 


y- 


Proof. We start with y e ynadm- As lm(A) is dense in y, it exists a sequence 
£ (Y in such that Axn —-—> y. This sequence cannot have any bounded 

subsequence: indeed, if such a subsequence would exist, there would be another 
subsequence, denoted x^ here, such that x^ weakly converges to an element x in 
At. The operator A being linear and strongly continuous, it is weakly continuous 
[38], hence Axm weakly converges to Ax. But by definition Axm strongly converges 
to y. By uniqueness of the limit, we have Ax = y, and y £ yadm, in contradiction 


with the initial assumption. Therefore, we have ||2:^11 a 


■A +00. 


Now, consider y £ yadm- The previous proposition implies the existence of a 


sequence ym £ ynadm such that yr, 


-A y. For a fixed m, we now know the 


existence of a sequence Xm,n £ A such that Ax 


n—>-oo 

m,n 7! ^ 2/^ 




and Wxri 


I n—^oo 

U -^ 


+ 00 . In particular, for any m £ N, there exists n(m) £ N such that Xm '■= Xm,n(m) 
verifies at the same time 

1 


\x>'m and \\Axm - U. 


m\\y < 


m 


It is then not difficult to verify that the sequence verifies the researched prop¬ 
erties. □ 


Remark 1. Actually, if y is not an admissible data, it is shown in the proof that 
any sequence such that Axn - y verifies lim ||a;„|j;f = -l-oo. 

3 ^ n—¥oo 

This proposition has for important consequences the fact that for any admissible 
data y, with corresponding solution cc, one can find an admissible data y, with 
corresponding solution x, such that y is arbitrarily close to y and x is arbitrarily 
far from x. 

We retrieve here the well-known fact that the problem of noisy data is crucial 
in data completion problems. Clearly, it is not sufficient to build a method that 
(approximately) reconstruct the solution of the data completion problem for any 
admissible data, it is also mandatory to propose a strategy for noisy data, as in 
practice data are always corrupted by some noise due to inaccurate measurements. 


3. Standard quasi-reversibility method. We define b a symmetric bilinear non¬ 
negative form on At, and denote by ||.||;, the induced seminorm on At. We suppose 
that it exists two strictly positive constants c, C such that 

c2||x||^<P2:||5, + ||2:||^<C2||x||^. 

Therefore, the symmetric bilinear form (., .)A,b, define by 

V(a;, xt) £ X, (x, x)A,b = (Ax, Ax)y + b(x, x), 

is a scalar product on X, and X endowed with this scalar product is a Hilbert space. 
We denote H-IU.b the corresponding norm, which is equivalent to the ||.||a' norm. 
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Obviously, there exists such a form b: it suffices to take the whole scalar product 
in A", 6(.,.) = 

Adapting the initial idea of Jacques-Louis Lions and Robert Lattes [19], the 
quasi-reversibility method applied to the abstract data completion problem defined 
above relies on the resolution of the following regularized problem 


Problem, for y G y and e > 0, find G X such that 

{AXe, Ax)y + e b{Xe, x) = {y, Ax)y, Vz € X. 

The quasi-reversibility equation is the Euler-Lagrange equation corresponding to 
the minimization over X of the energy \\Ax — f\\y + e||a;||j. In other words, it is a 
Tykhonov regularization of the data completion problem, e > 0 being the parameter 
of regularization and the penalization (semi)norm. Since its introduction in 
1963 by A.N. Tykhonov [37], this regularization has been widely studied and used to 
solve inverse problems (for a complete study on the topic, see [36] and the references 
therein). There are various methods to study such regularization method: e.g. 
singular value decomposition if A is compact (which is not the case in our data 
completion problems, see section 5) or spectral theory. Here we propose another 
approach to study the method, based on the variational formulation of the quasi¬ 
reversibility method, and on the differentiability of the approximated solution with 
respect to the parameter of regularization, the later being useful in the study of the 
iterated quasi-reversibility method. 

First of all, let us verify that the quasi-reversibility problem is well-posed. 


Proposition 3. For any y G y and e > 0, the quasi-reversibility problems admits 
a unique solution x^, with the following estimates: 


IIAxelly < ||/||y, \\Axe-: 


b< 


bell^ 


i(l, Ve) 


lly|| 


V- 


Proof, let us define the bilinear form 

ae{x, x) := {Ax, Ax)y -I- e h{x, x), Vx, x G X. 

It is obviously continuous. Furthermore, for all x G X, we have 

ae{x,x) > min(l,£:)||a;|b ,„ 

and therefore it is coercive. Finally, as \{y,Ax)y \ < \\A\\ \\y\\y \\x\\x < \\A\\ \\y\\y b|U.b, 
we obtain the existence and uniqueness of Xg by Lax-Milgram theorem. By defini¬ 
tion, we have 

\\Axe\\y < \\Axe\\y + e\\xe\\l = {Axe,y)y < \\Axe\\y\\y\\y => \\Axe\\y < \\y\\y. 
Furthermore, 

{Axe-y,Axe)y = -e\\xe\\l < 0 ^ \\Axe-y\\y < -{y,Axe-y)y < \\y\\Y\\Axe-y\\y, 
implying \\Axe - y\\y < bUy. Finally, we have 

\\Axe\\y -G e\\xe\\l = {Ax,,y)y < \\Ax,\\y\\y\\y < ^\\Axe\\l,-Ge\\xe\\l ||yb 
leading to min(l, ^/e)\\xs\\A,b < \J\\Axe\\y +£\\xe\\l < \\y\\Y- □ 


£ 


A' 


> 0 . 


Remark 2. In particular, we always have x, 


e—>-oo 
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Suppose there exists x G X such that Ax = y [i.e. y G yadm)- It is easily 
seen that x is never the solution of the quasi-reversibility problem, except in the 
special case y = 0 (which is always in yadm) for which x = 0 = x^- In other words, 
there is no £ > 0 such that the quasi-reversibility method reconstructs exactly the 
exact solution of the data completion problem. As seen in the following corollary, 
the solution of the quasi-reversibility problem is also never 0, except again in the 
special case y = 0. 

Corollary 1. The three following properties are equivalent: 

(i) y^O 

(a) 3 £ > 0 S.t. Xe ^ 0 
(Hi) V£ > Q, Xe ^ 0. 

Proof, obviously, (iii) implies (ii). Furthermore, as min(l, •y/£)||a;e||A,b < l|y||vj (ii) 
implies (i). 

Suppose it exists £ > 0 such that Xg = 0. For that particular e and for any 

- y 

X G X, we would have {y,Ax)y = {Axe,Ax)y +e{xe,x)b = 0. As Im{A) = y, 
this directly implies y = 0, so (i) implies (iii). □ 


Proposition 4. Lety G y, andxg the solution of the corresponding quasi-reversibility 
problem. Then Ax^ strongly converges to y (even if y is not an admissible data). 


Proof. As min(l, •\/£)||a;em_6 < ||y||y, we have, for any x G X 

^ . £—>-0 


e\b{xe,x)\ < 


min(l, 


llyl|2;||b 


> 0 . 


Let (£m)meN be a decreasing sequence of strictly positive real numbers such that 
limm->.oo £m = 0, and note Xm ■= As |jAa;m||y < ||y||y, it exists a subsequence 
(still denoted Xm) such that Axm weakly converges to y G y. But, for all x G X, 
we have 


(y - y, Ax)y yXATX (y _ Axm, Ax)y = emb{xm, x) q, 

- y 

that is y = y as lm(A) = y, and Axm weakly converges to y. As ||Aa;m||y < WuWy 
(proposition 3), Axm strongly converges to y. It is then not difficult to see that 
Axg strongly converges to y as £ goes to zero. □ 


We can now state the main theorem regarding the standard quasi-reversibility 
method: 


Theorem 3.1. Suppose y G yadm, and let Xg be the (necessarily unique) solution 
of the abstract data completion problem. Then Xg converges to Xg as e goes to zero, 
and we have the estimates IjAxe —yjjy < -yiUxsIlb, ||a::e||b < l|3:^s||6 o-nd H^e —Xs||b < 
ll*s||b- 

Suppose y G ynadm- Then lim ||xe||h = -l-oo. 

£— ^0 

The theorem remains valid when the ||.||b seminorm is replaced with the H-llA.b 
norm. 


Proof. Suppose hrst that y G ynadm- Then, as Xg is a sequence in X such that Axg 
converges to y (proposition 4), proposition 2 and remark 1 imply limE_j.o Ha^ellA = 
-hoc. As ||a;,,||^ j = WAxgW^ \\xg\\l > c^\\xg\\x, we have lim,,^o lla^ellb = +oo. 

Now, suppose it exists Xg such that Axg = y. Then, choosing x = Xg — Xg as test 
function in the quasi-reversibility problem, we obtain 

\\Axg - y\\y eb{xg, Xg - Xg) = 0, 


( 1 ) 
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which in turn implies b{xs,Xs - Xg) < 0 ^ ||a;e||b < ||xs||;, ^ Ha^elU.h < lla^slU.fe- 
Therefore, Xg is a bounded sequence in X, and up to a subsequence it weakly 
converges to x. As A is a linear continuous operator, and hence is weakly continuous, 
proposition 4 implies Ax = y, which implies a; = ccg as A is one-to-one. The 
uniqueness of the limit implies that the whole sequence weakly converges to Xg. 
Finally as HxeUa.p < lla^slU.pi the sequence strongly converges to Xg- 
Subtracting £b{xg,Xg — Xg) to equation 1, we obtain 

\\Axe - y\\y + ellxs - XeWl = -eb{Xg,Xe - Xg) ^ lliCs - XeWl < \b{Xg,Xe - a:s)| 

and by Cauchy-Schwarz inequality, ||a;e — XsUb < ||a;s||b- 
Finally, equation 1 implies 

- y\\y < elkellb \\Xe - XsUb < e||Xs||t 

which ends the proof. □ 

Next, we focus on the differentiability of the solution of the quasi-reversibility 
method with respect to e, a result that will be useful in the study of the iterated 
quasi-reversibility method. 

3.1. Differentiability of the quasi-reversibility solution with respect to e. 

It turns out that Xg, solution of the quasi-reversibility problem, depends smoothly 
on the parameter of regularization e. Indeed, let us define the map F : e > Q ^ Xg. 

Proposition 5. The map F is continuous. 

Proof. We choose e > 0 and h such that e — \h\ > 0. For any x G X, we have 

(Axg+h, Ax)y + {£ + h) b{xg+h, x) = {y, Ax)y 
{AXe,Ax)y +£ b{Xe,x) = {y,Ax)y. 

Subtracting the two equations, and choosing x = Xe,h '■= Xgj^h — x^, lead to 
II AXg^ji II ^ -f II lib — bi b(^Xg;.i-fi , Xg^fi). 

In conclusion, we have 

min(l,e) Hi^.blllb < h\\xg+h\\b\\xE^h\\b < {e + h)~'^/'^)\\y\\y\\xg:^h\\A,b 

which ends the proof. □ 

Remark 3. If the data completion problem admits a solution Xg, then F extends 
continuously to K+ by defining F(0) = Xg. 

Proposition 6. IFe have F G (7^(151+; A). For all e > 0, F'^e) = xi^\ unique 
element of X verifying 

{Ax^J'\Ax)y + eb{x^^^\x) = —b{xe,x), VxGX. (2) 

Furthermore, ||ai£ \\A,b < Diin(I, £“ 2 )||y||j;. 

Proof. By Lax-Milgram theorem, there exists a unique xi^'^ G X verifying 2, and it 
clearly verifies 

min(I,£)||4^)||i t, < llxellbll^^^lU.b < min(I,£-5) ||?/b||a;W|U.b- 
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It is a continuous function of e: indeed, for e > 0 and /i G M s.t. e — \h\ > 0, we 
have, for all a; G A", 

Ax)y + (e + /i) x) = -b{xe+h-, x) 

{Ax^^\ Ax)y + e b{xA\x) = —b{xe, x). 

Choosing x = x^^l := — x^e'^ G A” and subtracting the two equations lead to 

+epSllb = -hb{x^^lh,S:^el)-Kxe,h,x''^l)- 

Therefore, 

Pe^llU.b < min(l, (e + + min(l, (e + min(l, \\y\\y 

implying the continuity of the map 9 e i—>■ xi^'^ G X. Remains to be proved that 

F'(e) = x^e'^. For e > 0 and G M such that e —\h\> 0, we have 

{AXe+h, Ax)y + (e + /i) b{xe+h, x) = {y, Ax)y 
-(Aa;e, Ax)y - e b{xe, x) = -{y, Ax)y 
—h{AxA\Ax)y — he b{x^^\x) = hb{xe,x). 

Choosing x = x^^h '■= x^yh — Xg — hx^'^ and adding the three above relations lead 
to 

\\Axe,h\\y+£\\xe,h\\l = -hb{xs,h, Xe,h) => pe,b|U,b < h'^ C{h,e)\\y\\y, with C{h,e) > c> 0. 
The result follows. □ 


A simple induction leads then to the following theorem: 
Theorem 3.2. F G A"). For e > 0, for all m G N, 

dFF 


de^ 


-is) := 4 -) 


with defined recursively by 


(0)_ 

Xs .— X^^ 


Vto G N, is the only element of X verifying 

, Ax)y + e , x) = — (m + 1) b{xi^'’ , x), Vx G A”. 

In particular, Xe™^ verifies the following estimate: 




m ! 




min(l, e”i+i/2) 

If the data completion problem admits a solution Xg, it is not difficult to prove 
that 

|U,b < min(l, e™)“^ m! ||Xg|U,b- 
Finally, we have the following generalization of corollary 1: 


Corollary 2. the three following properties are equivalent: 

(i) y fit) 

(a) 3e > 0, 3m G N s.t. xi™^ yh 0 
(Hi) Ve > 0, Vm G N, xi™^ ^ 0. 
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Proof. Clearly (iii) implies (ii). Furthermore, as min(l, < ml Hylly, 

(ii) implies (i). 

Suppose it exists £ > 0 and m € N such that =0. If m = 0, then corollary 1 
implies y = 0. If m > 0, as = —m\\x^J^~^'^\\l, 

we obtain xi™ = 0, and by induction x^ = 0, implying again y = 0. Therefore 
(i) implies (iii). □ 


3.2. Monotonic convergence of the qnasi-reversibility method. In this sec¬ 
tion, 2 /^ 0 . Using the results on the derivatives of x^ with respect to e, it is easy to 
prove that if the data completion problem admits a solution Xg, then Xg converges 
monotonically to Xg when e goes to zero. This is of course not the only method to 
obtain such results (see for example [36], where spectral theory is used), but it has 
the advantage to be quite simple. 

The main result of this section is the following 


Theorem 3.3. Suppose the data completion problem admits a unique solution Xg. 
Then ||xe — Xs||A,b is strictly increasing with respect to e. 


We need to prove first the following two results, which are true whether or not 
the data completion problem admits a solution: 

Lemma 3.4. For all m S N, for all u G N, (— l)’”~'’" 6 (x^™\ x^"^) > 0. 

Proof. For m € N, let us define the axiom of induction: 

P(m) : Vn G {0, • • • , m} , (-I)™+” 6 (x('"\ x^")) > 0. 

Obviously, P(0) is true, as j/ ^ 0 ^ x^ 0. 

Suppose P{M) is true for some M G N. Let fc be in {0, • • • ,AI + 1}. 

• if fc = M -f I, then > 0 (as y 7 ^ 0) 

• if fc = M, then, by definition of 


(-l)2^+i6(x(^+i),xW) = -&(x("^+i),xW) = 


\\Ax. 


(M-H)||2 


-I- e x 




> 0 


M + \ 

• if fc < M, then, using successively the definition of and we 

obtain 


5(x, 


(M-H)^3,(/c)) ^ 


-1 
fc -I- 1 




Ak+l)y 


As fc -I- I G {0,..., M}, P{M) implies 

(-l)“+''+^6(x(^),xf+i)) > 0 ^ (-l)“+'=+^6(x(^+^\xf)) > 0. 


□ 


Proposition 7. The quantity ||Ax£ — y||y is a strictly increasing function of e. 
Proof. Defining y : e G i-G- ^HAxg — y\\y, we have 

g'ie) = {Axg - y, x(^^)y = -£b{xg, x^^^) > 0 . 

□ 
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Proof of theorem 3.3. define g := e G !->■ ^||a;£ — Xs||b. We have g'{e) = b{xe — 
Xs.x^e'^) and 5 "(e) = \\x^e\\l+ h{xe - Xs,xf'’). Therefore 
eg”{e) = e||x(^)||^ + £b{xe - Xs,xP) 

= e||x^^^||j — (A(xe — Xa), Ax^^^)y — 26(xe — Xs,x^^^) ( definition of x^^^) 

= e||x^^^||^ + £b{xe,x^^^) — 2b{xe — Xs,x[^^) ( definition of Xg and Axs = y) 
= s\\x^P\\l + £b{xe,xi^'>) - 2g'(e). 

So g verifies the following ODE: £g"{£) + 2g'(e) = £||xi^^||^ + e&(xe,xi^^), that is 

^{e^g'ie)) = £^\\xP\\l + £H{xe,x^^^) > 0. 

Therefore, e^g'ie) is a strictly increasing function. As \\xa — Xs||;, < ||xs||{, and 

—rlla^sllh, we have 
min(l, e) 

|e^g'(£)| = |£^6(Xe - Xs,x(^^)| < £^||Xe - Xs||b||Xe||h < £||Xs||b -S> 0, 

£—>■0 

which leads to e'^g'is) > 0 ^ g'{£) > 0, which ends the demonstration, as ||xe — 
2 ^s||A,b = ^||Ax£ y\\y T ll^e ®s|lb' ^ 


4. Iterated quasi-reversibility. As seen in the previous section, the quasi-reversibility 
method can be viewed as a Tykhonov regularization of our abstract data completion 
problem. Therefore, it seems natural to study a well-known extension of such reg¬ 
ularization, namely the iterated Tykhonov regularization method, to our problem: 
we then obtain the iterated quasi-reversibility method. 

The iterated quasi-reversibility method consists in solving iteratively quasi-reversibility 
problems, each one depending on the solution of the previous one. More precisely, 
we define a sequence of quasi-reversibility solutions by induction : X~^ = 0^ and 
for all M € N, is the unique element of X verifying 

{AX^, Ax)y -\- eb{X^, x) = {y, Ax)y -G £b{X^~^,x), Vx G X. 

It is not difficult to verify that the sequence is well-defined. In particular, it is clear 
that A° = Xg, solution of the quasi-reversibility problem. 

Our study of the iterated quasi-reversibility method is based on the follow¬ 
ing result, which highlighted the link between the solutions of the iterated quasi¬ 
reversibility method (Xj'^)MG{-i}uN and the derivatives of Xg with respect to the 
parameter of regularization e: 


Theorem 4.1. For all £ > 0, for all M G { — 1} U N, we have 


M 


xf = y {-IT—, 

^ m ! 


Am) 


m—0 


Proof. Denote X^ ■= {—^T—^ ~ empty, therefore 

m—0 ^ ' 

we have X~^ = Ox = For M = 0, we also have X^ = x^ = A°. Finally, for 
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M >1 and l<m<M+l, in virtue of the definition of we have for all a; € 

Summing for m = 1 to M + 1, and adding the equation verified by = x^, we 
obtain 

{AX^+^, Ax)y + £b{x^+'^ ,x) = {f,Ax)y + £b{X^,x). 

A straightforward induction ends the proof. □ 

From now on, we suppose y ^ O', if not, we have X^ = 0 for all e and M. 

4.1. Some estimates on X^ and AX^. We start with estimates on the M-th 
iterated quasi-reversibility solution, valid for any data y, admissible or not. In other 
words, these estimates are valid whether or not the data completion problem has a 
solution. 


Proposition 8. For all e > 0, for all M G N, we have 

(a) ||Af-i|l,<||Af||, 

(b) IIAAf b < \\y\\y 

(c) \\AX^-y\\y<\\y\\y 

(d) \\AX^-y\\y<\\AX^-^-y\\y. 


Proof. We start with estimate (a): as y 7 ^ 0, we have 0 = \\X^ ^||t, < ||a;e||f, = ||Alg Hh. 


Furthermore, for M gN, ||Aj^||^ = EE(- Therefore, 

from lemma 3.4 we obtain 


k=0 m =0 


I A. 


M+li 


- IIA, 


Ml 


M-l-1 

= 2 ^ (_l)M+l+m 

m—0 




{M -I-1)! to! 


b{x 


(M+l) m 


)> 0 . 


Regarding estimates (b) and (c), we note that they hold for M = 0. Furthermore, 
\\AX^+% = (y,AAf+ib+£6(Af,Af+i)-e||Af+i^ 

Estimate (a) implies b{X^~^^,X^) < and 

\\AX^+^fy<{y,AX^+%. 

Cauchy-Schwarz inequality implies then the estimate (b). Furthermore, 
\\AX^+^-yfy = iAX^+^-y,AX^+^)y-{AX^+^-y,y)y < -(AX^+^-y,y)y 
which leads to estimate (c). 

Finally, the case M = 0 of estimate (d) correspond to estimate (c) with same M. 
For M G N, we note that 

M-l-l ^ 

IIAAf+1 - yfy = II ^ (-1)™^a4-) - yfy 

m—0 

o p A/+1 

= PXf - yWl + ylxf - y)y. 
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Therefore, it is sufficient to determine the sign of — y)y. 

By definition, we have 

{Axi^+^\Ax,-y)y = -e b{xi^+^\x,) 
and for all m G { 1 ,..., M}, 

^m +1 




Nm+l 


+ (- 1 )™+! 




(to — 1)! 




Summing these equations for to = 0 to M, we obtain 


iAxi^+^\AX^ - y)y = 


pAf+1 

^ ^ (M+1) 


M\ 


b{x\ 




In conclusion, {—l)^'^^{Axi^^^\ AXf^ — y)y = ^j^b{xi^'^^\xi^'^) < 0 by 
lemma 3.4. The result follows. □ 

Proposition 9. For all e > 0, for all M G { — 1} U N, ||Xj^||f, < 


2(M+1) 


bll: 


Proof. Proposition 9 is obviously true for M = —1. Let M G N. We consider first 
the inductive sequence: 

xi > 0 

VM G N, xm +1 > 0 and — xm+iXm — xf = 0. 

Note that the sequence is well defined as pM(a;) '■= x^ — xm x — x\ verifies pm( 0) < 0 
and therefore has a unique strictly positive root. 

We prove by induction that xm < V‘2 Mxi. It obviously holds for M = 1. 
Suppose that xm < \/2Mxi for some M gN. Then 

Pm (^\/2{M + l)a;i^ = 2(M + l)x\ - a/2(M + 1)xiXm - x\ 

> (2M + l)xl - 2VM + iVMxI 

= (yM + 1 - Vm^xI > 0 

and therefore xm+i < \J2{M + \)xi. 

Now, we specify the sequence, defining xi := ^p||?/|jy, and prove by induction 

y/e 

that ||Xj'^||b < xm+ 1 - Suppose it holds for some M G N. Note that by definition 
of we have 

||GlXf+i|| 2 , +£||xf+i ||2 = {y,AX^+^)yFeb{X^,X^+^) 

<\\y\\y\\AX^+^\\y+e\\X^UX^+% 

< \\y\\y + exM+i\\X^‘"^^\\b, 

(we used estimate (b) of proposition 8 here) which in particular implies 
\\x^""^\\l - XM+i\\Xf+^\\b -xl<0. 

The definition of the sequence (xm)mgn implies directly ||X^+^||{, < xm+ 2 - As the 
result is true for M — 0, the proposition follows. □ 


Remark 4. Actually, for M G N, the inequality in proposition 9 is strict. 
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4.2. The case of exact data. From now on, we suppose that y S Yadm, and denote 
Xs the solution of the abstract data completion problem. We define := — Xs, 

the discrepancy between the M-th iterated QR solution, and the exact solution. 
Note that by definition, = —Xs and for all M € N, 

{AR^, Ax)y + £ b{R^, x) = e b{R^, x), \/xGX. (3) 

We aim to prove the following theorem: 


Theorem 4.2. For all e > 0, R^ 0. In other words, for any £ > 0, 

converges to Xg as M goes to infinity. 

As X^^ = it means that the sum converges as M goes to 

infinity. In other words, it means that if it exists Xg S A solution of Axg = y, then 


Xs 


Hi-i) 

m—0 


m ! 


.(m) 

e > 


hence the solution of the data completion problem can be seen as a series of deriva¬ 
tives of the quasi-reversibility solution w.r.t. the parameter e. 

Let £ > 0 be fixed. We start with the following estimates 


Proposition 10. For all £ > 0, for all M G { — 1} U N, 

- \\R^+%<\\R^U 

- I!^fll.<l|xg||,. 

As a consequence, we have ||i?^||b < ||xs||6 for all M G {—1} UN and all £ > 0. 


Proof. Choosing x = R^ in (3), we obtain 


IIARf ||^+£||i?f 11^ = £6(Rf-\Rf) ^ WR^Wl < 6(i?f < ||Rf ||,||Rf "i 

hence the first estimate is valid. Note that in particular, as |ji?^||b < llRellb = 


b, 


II Xe — Xs||b, we have ||i?^||b — 0 for any M gN. 

Let us now focus on the second estimate, which is directly true for M = — 1. 
For M e N, let us define g := e G Kj i-)- |||Arj^||^. As by definition, X^ = 

have 



M 


Hui) 

m—0 


(m+1) 


M 


Eui) 


£m-l 

(m — 1)! 




~1V1 

^ > M\ ^ 


Therefore, we have 


g'{e) = b{-Xf,Xf) = 


M 

m—0 


M-\-m 


M\\ 


-6(xi^+i),x(" 


< 0 , 


implying in particular that ||Xj'^||b < lim ||Ar,^||b = ||xs||&, as |li?£^||h XAffi. q. □ 

ry—^0 ' 


Proposition 11. The series ||Aii^||y converges, therefore 

M 


lim 

Af —>-oo 


||Ai?fb = 0. 
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Proof. For any M € N, we have 

\\AR^^+Yy + = eb{R^ 

which leads to 

\\AR^+Yy < eb{R^,R^+^) < e||i?f ||,||i?f < e\\R^\\l = e6(i?f-\i?f )-Pi?f ||^. 
Therefore, we obtain 

\\AR^+Yy + lly < eb{R^-\R^) 

and by an immediate induction 

M 

iT\\y<4RX<4^s\\l 


E 11^^” 

m—1 

Therefore , the series E] P^e"lly converges. The property follows. 


□ 


Theorem 4.2 follows from the previous proposition: indeed, let : N —>■ N be a 
strictly increasing map, and define R^ := R!§^^\ As 

||i?f Hi,, = WARffy + WRfWl < lli/lli + WxsWl 

we have that {R^)MGn is a bounded sequence in X. Consequently, there exists 
1? : N —>■ N, a strictly increasing map such that R^ := weakly converges to 

R°° in X. As A is linear and continuous, we directly obtain from proposition 11 
that AR°° = Ov, which implies R°° = Oa- as A is one-to-one. 

We hence have obtained that R^ Ox, or equivalently ^ Xg- But 

Af->-oo M^oo 

we know from propositions 8 and 10 that < ||xs||a,&, implying that 

strongly converges to Xg, and it is then not difficult to show that the whole sequence 
Xy strongly converges to Xg as M goes to infinity. 


4.3. The case of noisy data. In this section, we suppose that our exact data, 
denoted yex G y, for which the data completion problem admits a unique solution 
Xg G A, is corrupted by some noise. The obtained perturbed data, denoted G y, 
is supposed to verify ||?/'^ — yex\\y < in other words, we know the amplitude of 
noise on the data. On the other hand, there might or might not be x G A such that 
Ax = y^: we don’t know if y^ is an admissible solution or not. 

From now on, for any y G y, we will denote X^{y) the M-th iterated quasi¬ 
reversibility solution with y as data. Our main objective in this section is to propose 
an admissible strategy to choose M as a function of <5, the amplitude of noise, to 
ensure that, when 6 goes to zero, tends to the exact solution Xg. As pointed 

out in proposition 2 and remark 1, this is a crucial point in the study of data 
completion problems. 

A first important remark is the following: AX^{y) always converges to y, re¬ 
gardless of the admissibility of y as data for the data completion problem. 

Proposition 12. For any y G y, for any e > 0, AX^{y) — - — y. 

M—>.oo 

Proof. As jVadm IS dense in y, for any r; > 0, it exists yrj G yadm such that \\yjj — 
y\\y < 5- Proposition 8 (b) implies ||AXy(2/) - AXy(y^)||y < \\y - y-,^\\y < 
Finally, as y^ G yadm, there exists Mn > 0 such that for any M > 
||AXy(j/^) — y.qWy < y The result follows. □ 
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Next proposition defines the admissible choices of M to ensure the desired con¬ 
vergence: 


Proposition 13. For any e > 0, for any choice of M := M{6) such that M{5) 

-boo and 5^JM{5) - > 0, we have > x^. 

5—>0 (5—>0 


-S> 


Proof, obviously, we have, for any £ > 0 and M G N 

l!^f (/) - x,\\A,b < ll^f (/) - X^iy)\\A,b + ||Xf (y) - Xs\\A,b- 

If M := M{5) verifies lim5_>oM(^) = -boo, theorem 4.2 implies directly that 

\\X^'^^\y) — cCsIU 6-S' 0. Furthermore, propositions 8 and 10 imply 

’ 5—>0 


The result follows. 


1 


2 (M-b 1) 


)ll/-y|lv = (i + 


2 (M-b 1) 




□ 


Proposition 13 defines the admissible strategies to choose M depending on 6. 
An admissible choice could be M{5) := for example. But such a choice, if 

it guarantees the convergence of the method, does not correspond to any precise 
objective. We therefore focus on another method to choose M{5). 

Letr > 1. Forafixede > 0, we define := {M G N, \\AX^ {y^) — y^\\y < rJ}. 
Proposition 12 implies that is non-empty, and we define M (5) as the minimum 
element of M^: M{5) := min{M G M((5)}. 

M{5) is chosen accordingly to the Morozov discrepancy principle: it is the first 
M G N such that the distance between AX^ and y^ is (approximately) equal to 
the distance between Axg = y and — y^Wy ~ ll^la^s — y^\\y- This method 

to choose M depending on S has two interesting characteristics: 


1 . with this choice, one does the minimum number of iterations of the iterated 

quasi-reversibility method required to obtain an error in the residual || — 

of same order of the error on the data. 

2. such choice is admissible, in the sense of proposition 13. 

We now prove that M{S) is an admissible choice. 


Proposition 14. M{S) — -boo. 


Proof. Suppose it is not the case. Then there exists a sequence of strictly positive 
real numbers Sn and a positive constant C such that 0 and M{5n) < C. 

It implies the existence of a subsequence (still denoted Sn) and Moo G N such that 
Sn 0 and M(Sn) Moo. In particular, it exists N G N such that for all 

n>N, M(Sn) = Moo- 

For n> N, the definition of M{S) implies 


A 0. 


||AAf»(/")-ye.b < \\AX^-{yi)-yHy + \\y^--ye.\\y < {r + l)Sn - 

Consequently, using proposition 8 we have 

||AAf-(ye.) - ye4y < \\AXf-{y,4 - AXf-{44\\y + ||AAf~(/") - y,4y 


< hex - y^'Wy + (r + hh = (r + 2)Sn 


-^0, 


n—>-oo 
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that is AX^°° iVex) = Vex- If we directly obtain y^x = 0, in contradiction 

with the hypothesis. If M^o = 0, we obtain = Xs, which again (corollary 1) 
implies ys = 0 Finally, if M^o > 0, as for all a; € d:”, 

{AX^°°{yex),Ax)y + eb{X^°°{yex),x) = eb{X^°°~^{y^x),x) + {yex,Ax)y, 

we obtain X^’^{yex) = X^^°°~^{yex), or equivalently xi^°°\yex) = 0, which again 
implies yex = 0 by corollary 2. We obtain once again a contradiction, which ends 
the proof. □ 

Proposition 15. lim 5^JM[5) = 0. 
i5—>-0 

Proof, by definition, we have \\AX^^^'^ ^iv^) ~ V^Wy > Therefore 

WAX^m-i^y^x) - yexWy = pxf W-I(yex) - - 

> pxf - y^Wy - pxf - 

> (r — 1)5 

as proposition 8, estimate (c) implies 

Pxf - AXf (')-!(/) - {y,x - /)\\y < hex - y^\\y < 6. 

We hence have obtained: 

{M{S)-l)Shr-ir < (M(5)-l)pxf = {M{S)-l)\\AR^^^^-\y,x)\\ 

As \\AR^+'^{yex)\\y < \\ARf{yex)\\y and \\AR^{yex)\\y converges, m\\AR^{yex)\\y 
tends to zero as m goes to infinity. Therefore, {M{5) — l)||Ai?^^‘^^~^(yea;)||y goes 
to zero as 5 tends to zero, implying that 

lim(M(5) — l)5^(r — 1)^ = 0. 

5->-0 

The result follows. □ 

5. Quasi-reversibility methods for data completion problems for the Pois¬ 
son’s equation and the heat equation. We will now go back to the data com¬ 
pletion problems described in the introduction, and verify that they correspond to 
the abstract setting introduced in section 2. 

5.1. Poisson’s equation. As mentioned in the introduction, the data completion 
problem for Poisson’s equation is: for {f,gD,9N) S x T^(r) x T^(r), find 

u G H^{n) s.t. 

{ —V • crVu = / in n 

u = go on r 

aXu ■ V = gv on T. 

We could directly use this formulation of the problem to obtain a quasi-reversibility 
regularization. However, if we do so, we obtain a fourth-order variational problem, 
which is rather difficult to discretize as we would need or non-conforming finite 
elements which are seldom available in numerical solvers. Therefore, we first modify 
the problem by introducing the flux p := cr-Vu as an additional unknown, following 
the idea introduced in [31]. It verifies —V • p = —V • crVu = / G L^(H) and 
p . 1 / = ctVu ■ V = gN & A^(r), hence 

p G #div(fl) := {q G L'^{nf, V • q G q • z/ G T^(r)} . 


yex\\y 

{.Vex ~ y ) lly 
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endowed with the scalar product 

(P, := ^ (p • q + (V • p) (V • q)) dx + J^{p • u) (q • u) dS 

is an Hilbert space [32] . We modify the data completion problem the following way: 

Problem. For f, go and gjq in respectively L?{T) and L‘^{T), find (u,p) G 

F[^(fl) X such that 

{ aVu = p m H 
—V • p = / in H 

u = go on r 
p-iy = g^ onT. 

Obviously, this is exactly the same problem as previously. However, this small 
modification will lead to a second-order variational quasi-reversibility regularization 
in the product space x Hd\v, easily discretized using standard finite-elements. 
To fit in our abstract setting, we introduce the operator 

A :(u,p) G A = H\n)x 

(crVn - p, -V • p, U|r, p • v\t) G 3^ = x L^{n) x ^^(r) x ^^(r). 

The spaces X and y, endowed respectively with the scalar products 

((M,p),('y,q))^ := {u,v)m + (p>q)H,i„ 

and 

(^{F,f,g,h),(F,f,g,h)^^:=J^(F-F + ff^dx + j (^gg + hh^ ds 

are obviously Hilbert spaces, and the data completion problems can be rewritten: 
find (m,p) G X such that A{u,p) = {0, f,gD,gN) G y. 

Proposition 16. A is linear, continuous, one-to-one. It is not onto but has dense 
range. Additionally, A is not a compact operator. 


Proof. Clearly, A is linear continuous. As the data completion problem for Poisson’s 
equation is known to admits at most a solution, but may have no solution, A is one- 

- y 

to-one but not onto. Let us prove that lm(A) = y. Let (F, f,g, h) G y such that 
(^A{u,p),{F,f,g,h)^^ = 0, y{u,p)Gy, 

that is 

(crVu—p)-F—(V-p) /) dx-l- / (ug+{p-i') h\ ds = 0, y{u,p) G H^{fl)xHdiv{^). 


Choosing u = (p G C'^{Fi) and p = 0, we obtain —V • cr^F = 0, and in particular 
cr'^F G Hd„{Fl). Choosing w = 0 and p = G we obtain V/ = F. 

Therefore, / G H^{Fl), and verifies —V • cr'^V f = 0. Hence, taking u G Il^(fl) and 
p = 0, and using the Green formula, we obtain 

(tT^V/-u,M)= j f-Wudx = I a\tu ■ f dx = — I ugds, Vw G 

JQ. Ja Jr 
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implying that cr^'^f ■ i/ = —g on T and 0 on Fc- Taking u = 0 and p G iJ(jjv(n), 
and using the divergence theorem, we obtain 

(P • i",/> = ^ (p • V/ + V • p/) da; = J^{p ■ v)hds, 

and therefore f = h on T and 0 on Fc- We have obtain that / G verihes 

—V • cr^V/ = 0 in O and / = cr^'^f ■ v = Q on Fc: by uniqueness of the solution 
of the elliptic data completion problem, necessarily / = 0, which implies directly 
F = 0 and g = h = 0. 

Finally, let us prove that A is not a compact operator. Consider e„ an Flilbert 
basis of and in id^(O) verifying / Undx = 0, —V • aVun = e„ in 

Jo. 

VL and aVun • = 0 on dfl. It is not difficult to show that Un exists and is 

unique. Furthermore, HunllErqn) < and in particular ||M„||^ 2 (r) < C(p.,a). 

Defining p„ := a'^Un G we obtain ||Pn||//j,^ < hence (wn,Pn) is 

a bounded sequence in X. But A{umPn) = (0, e„,|r, 0) does not admit any 
convergent subsequence. □ 


To define our quasi-reversibility approach to this data completion problem, we 
choose 6 (.,.) such that the corresponding norm ||.||A,b is equivalent to the norm 
||.||a'- Of course, we could choose the whole T’-scalar product. But it might be 
interesting to use another form, to soften the regularization: here we define 


&((u, p), (v, q)'j := Jjyu • Vu -F p • q) da; 

which is a symmetric bilinear non-negative form in X (but obviously not a scalar 
product). Using Poincare inequality, it is easy to obtain the existence of c, C > 0 
such that 

c||i;,q|U < lk,q|U,b < C'lk,q||x- 

As for any (u,q) G X, ||u,q||6 < ||w, qlla-, for a fixed e the regularization term in 
the quasi-reversibility method is smaller. 

Applying the abstract setting to this problem, we obtain the following quasi¬ 
reversibility regularization: for e > 0, find (ue,Pe) G id^(D) x il^div(n) such that for 
all (u,q) G H^{n) X Iddiv(U), 



;) • (crVu - q) dx -F J (V • Pe) (V • q) dx -F J (u^ "c -F (Pe • v) (q • 


ds 


+e 


(Vme • Vx -F Pe ■ q) dx 


/ (V • q)dx -F 


gov + gN {q- >^)) 


ds. 


JQ an ar ' ^ 

This problem always admits an unique solution (ue,Pe). We know from our study 
that if the data completion problem for the Poisson’s equation admits a solution 
(ms,Ps), then (Me,Pe) converges monotonically to (ms,Ps) as e goes to zero, with 
the estimate 



PeW 


2 

L^{Q)d- 






< Ve||Ms,Ps||&. 

If not, we know that ||(M£,Pe)||{,--Foo. 

6 —^0 

The quasi-reversibility method we obtain in this study is close to the one proposed 
in [31] to stabilize the data completion problem, which was: find {ug, Pe) G id^(U) x 
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-ffdiv(n): Us = qd and • u = on F, such that for all (u,q) S H^{n) x 
u = q • u = 0 on r, 



Pe)-(crVi; - q) dec -I- / (V • p^) (V • q) dx 
Jn 

+ e ((ue,Pe),(u,q)) =-/ /(V-q)da;. 

^ ^ ^ Jq 


The two differences are first the use of 6 (.,.) instead of {■,-)x in the regularization 
term, and secondly in the way the boundary condition are included in the problem. 
In the formulation proposed in [31], they are strongly imposed, which presents two 
main issues: one is theoretical, as the regularized problem might not have solution if 
qd is in T^(r), and not in iF^/^(r), as in that case there is no u € such that 

V = go on r, and therefore Ug cannot exist. The second one is practical: it is not a 
good idea to strongly impose data which might extremely noisy, as in that case the 
noise is somehow imposed to the solution. In the quasi-reversibility regularization 
obtain in the present paper, the boundary conditions are weakly imposed, which 
solves both of the problems: the regularized problem always admits a solution, 
even in the case where go is not the trace on F of a function, and the noise is 
regularized directly by the formulation, leading to a stabler formulation. 

Finally, the abstract iterated quasi-reversibility method applied to the elliptic 
data completion problem is: for e > 0, (n“^,p“^) = (0,0) and for all M G N, 
(wf ,Pf) e H\ri) X Hdiy{n) verifies 

[ {a\/uf - pf) • (aVn - q) dx + [ (V • pf ) (V • q) dx 

Jq J q 


+ ^ (wf + (pf ds + eb({uf, pf), (n, q)) = 

- ^ / (V • q)da;-h J 'c + ffiv (q • u)) ds-h £ Pf "f, (u, q)^ 

and we directly know that and p^ converge to u and ctVu when M goes to 
infinity. In the case where noisy data f^, g^ and gf^ are available, such that 

IlflD ~ 5i)lli2(r) + IldM ~ 5Alli2(r) + 11/ ~ /^Ili2(n) ^ 

in accordance with the result of section 4.3, we stop the iterations the first time 
that \\A{u^, p^) — ( 0 , /^, g^jj, < rS, with r > 1 close to 1 . 

Remark 5. Actually, in the following numerical results, we use r = 1. 

5.2. Heat equation. As for the Poisson problem, we modify the data completion 
problem defined in the introduction, introducing the flux p := Vm as an additional 
unknown: 

Problem, find (it,p) m L^(0, T; 7F^(f2)) n iF^(0, T; L^(n)) x L^(0, T;/fdiv(f^)) such 
that 

( dtu-V Tp = f in Q 

\ Vu = p in Q 

I u = go on (0,T) x F 

[ p v = ort (0,r) X F 
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Again, we define 

A:{u,p) G A := L'^{0,T;H\n))nH\0,T;L‘^{il)) x L‘^{0,T; 

^ (Vu -p,dtU-\/- p, u\T, P • J^ir) G = A"(0, T; L^inf) x L'^{0, T; L^{n)) 
X L‘^{0,T;L^{T)) x L^{0,T; L^{T)). 

Here, the spaces X and y are endowed with their natural scalar products, re¬ 
spectively 

^(m, p), (u, q)^ '■= f f (dtudtv + S/u ■ S/v + uv+C^ ■ p) {S/■ q) + p ■ c^dx dt 

^ J Q do 

+ f ({p ■ v) {q-v) dS dt 
Jo Jr 


and 


((Fi, /i,5i, /li),(F 2 , / 2 ,52, h 2 ))^ '=11^ (^1 • F 2 + /l f 2 )dxdt 

+ J J (^9192 +hi h2^dSdt. 

It is then not difficult to verify the 

Proposition 17. A is a linear continuous. It is one-to-one but not onto, and has 
dense image. Furthermore, it is not a compact operator. 

Proof. We will just prove prove that A has dense range, as it is not difficult to be 
convinced that A is non compact, and the rest of the proposition follows directly 
from the definition of A, X and y, and the ill-posedness of the corresponding data 
completion problem. 

Let F G L2(o,T;L2(o)<i), / g g G T; L2(r)) and h G 

L^(0, T; L^(r)) be such that for all v G L^(0, T; iL^(H)) n iL^(0, T; L^(H)) and all 
q G L^(0,r;H^div(^^)), (^(v.q), {F,f,g,h)^^ = 0, that is 

f f - q)-F + {dtv-S/■ q) f'^dx + f f (v g + {q-n)h^ds = 0. 

First of all, choosing q = T G C'“(H)'^, we obtain F = V/, and therefore / G 
L^(0, T; So we have 

[ [ (i'^v - q)-S/f + {dtv■ q) f'^dx + f f (v g + {q-n)h^ds = 0. ( 4 ) 

For all q G L^(0, T; ifdiv(f^)), for almost all t G (0,r), we have 

^ (q • V/ + (V • q)/) da; = (q • n, f)H-^/Hdn),H^/Hdn) 
which leads by integration in time to 

/ (q ■ = ~ / / {q-v)hdsdt 

Jo Jo Jr 

that is / = —h on (0, T) x F and / = 0 on (0, T) x Lc- 
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Now, taking v = ip € C^{n) in (4) leads to dtf + A/ = 0 in (0,r) x ft. We 
see that V := (/, V/) € -ffdiv(Q), and we can apply the divergence theorem: for all 
v€ =L2(o,r;7Ji(f7))niJi(0,T;L2(o)), we have 


J • V)?; -f V • Vt,xV^dxdt = (V • v,v) 

that is, taking any v S such that w(0,a;) = v{T,x) = 0, 

/ / + Af)v+f dtv+V f-Vvjdxdt = ((V/-U, i;)^-i/ 2 (aQ) ^/i/ 2 (aQ))^-i/ 2 (o,T),Hi/ 2 (o,T) 

Jo Jn -V-" ^ 


=0 


leading to V/ ■ i/ = —g on (0,T) x T and V/ • u = 0 on (0,T) x Fc- Therefore, / 
verifies dtf + Af = 0 in (0, T) x Q, f — V/ • u = 0 on (0, T) x Fc, hence / = 0 in 
(0, T) X F, leading to F = 0 and g = h = 0. □ 


Similarly as for the previous regularization, we introduce the symmetric bilinear 
non-negative form 

b^{u,p),{v,q)^ ■= ( f {dtudtv + \^u ■'Vv + p ■ c^dx dt. 

It is easy to check that the bilinear form {A{u, p), A{v, q))y + 6((u, p), (n, q)) is a 
scalar product on A”, and that it exists two constants c, C > 0 such that 

c||m,p||a < ||u,p|U,b < C\\u,p\\x. 

The quasi-reversibility regularization we consider is therefore: for e > 0, find 
(ME,Pe) G X such that for all (r',q) G X, we have 

J J (^{dtUs - V • Pe){dtV - V • q) + (Vue - Pe) • (Vu - q)'^dxdt 

+ f f (^UsV + {pe-i'){q-v)^dSdt + e f f [dtUe dtv + Ue ■v + Pe ■ dx dt 

J 0 'J T' 

T 


= j J f {dtv-V ■ q)dxdt + J J (^gDV + gNiq-i')^dSdt. 

According to our present study, this problem always admits a unique solution 
(ME,Pe) that converges to (it, Vm) when e goes to zero. The corresponding iter¬ 
ated quasi-reversibility is: (u“^, p~^) = (0, 0) and for all M G N, {u^, p^) is such 
that for all {v, q), 

^ ^ ((5iuf - V • p^){dtv - V • q) + (Vuf - pf) • (Vt; - q))dxtit 


+ f f v + {p^ ■iy){q-i')jdS dt+e f f (dtu^ dtV+Wu^-Wv+p^^-q^dx dt 

= j J f {dtv — ■ q) dx dt + J J ^go v + gN{q->^)'^dS dt 

+e [ [ {dtu^-^dtv + \7u^-^-S/v + pf-^ ■q)dxdt. 

Jo Ja ' 


6. Numerical results. 
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6.1. Elliptic equation. We consider a domain C with exterior boundary F 
and interior boundary Fc defined by 


dr := 

with 



cos(t) 

sin(t) 


t G [0, 2 7r] 


dr,:= 



cos(t) 

sin(t) 


t e [0, 2 tt] 


r{t) = 1 + 0.1cos(2t) — 0.05sin(3t), rdt) = 0.5 — 0.02cos(t) + O.lsin(t). 


We consider the problem of reconstructing a Robin coefficient rj on dTc from the 
knowledge of a noisy Cauchy data on i9F. Mathematically, we want to 

reconstruct a function u G and a function rj G C^(Fc) such that 


—Am 

= 0 

in n 

u 

= 9i 

on F 

di,u 

— 9n 

on F 

di,u + gu 

= 0 

on F, 


The Cauchy data igjj,gff) G T^(F) x T^(F) is supposed to correspond to an exact 
data {gD,gN) corrupted by some noise of amplitude 6 : 

llffi) ~ goWt^^T) + \\9n ~ 5iv|li2(r) ^ 

Our strategy to reconstruct g is therefore to compute and approxima¬ 

tions of u and Vrt with the prescribed noisy Cauchy data on F and no data at all on 
Fc, using the iterated quasi-reversibility method for the Poisson problem. Then, we 


obtain an approximation of g on Fc by simply taking the ratio % = — 


m(6) 

Ps ’I 


In our experiments, g = 0.5 -I- 0.3 sin(2 {9 — 57r/4)), 9 being the polar angle of a 
point X on Fc, and gjq = 1. The corresponding Dirichlet data is obtained by solving 
the direct problem — Au = 0 in fl, = 0.2 on F and d^u -I- ? 7 m = 0 on Fc using a 
finite-element method, and defining go := M|p. 



Figure 1. Exact solution. 

Then we corrupt the Dirichlet data go pointwise with a normal noise having zero 
mean and variance one, to obtained the corrupted Dirichlet data g^. The noise is 
scaled so that 

hi, - goWoo = ahoWoo 

that is the relative amplitude of noise in L°°-norm is a. In the experiments, we have 
chosen a = 1%, 2% and 5%. The exact Neumann data is used {i.e. gf^ = g^), as in 
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practical situations it is the imposed data (the net current), whereas the gu is the 
measured data (the corresponding voltages). Therefore is known quite precisely 
compared to go- We then compute the corresponding amplitude of noise for the 
norm S = S{a, ||gD||oo), which defined our stopping criterion for the iteration of 
the method. 

The iterated quasi-reversibility problem is then solved using a conforming finite- 
element method using P 2 Lagrange finite elements for and RTi Raviart-Thomas 
finite elements for [39] . The study of convergence of the finite-element approx¬ 
imation of the quasi-reversibility approximation toward the continuous solution is 
just a slight adaptation of section 4.4 in [31], as the formulations are quite similar, 
and therefore is omitted in the present study. To avoid an inverse crime, the direct 
and inverse problems are solved on different meshes. 



Figure 2. The direct and inverse meshes used in the numerical simulation. 


According to our study, the choice of e is completely arbitrary in the iterated 
quasi-reversibility method. Therefore, we have chosen e = 1 in the experiments, as 
it leads to a good conditioning of the finite-element matrices. 

First of all, we present in figure 3 the evolution of the residual 

^l|v • pf 11^2(0) + llVwf - pf Ili2(n) + \\Ue - gUhiD + IIPe • ^ - 9%\\lht) 

until the stopping criterion is reached. As expected theoretically, the greater is the 
noise, the smaller is M{5). 

Now we present the reconstruction results: in figure 4, the exact solution is 
compared to the reconstructed one in the whole domain of study. In figure 5, we 
focus on the boundary F: we compare the exact data gD^ the noisy one used 
in the iterated quasi-reversibility method, and finally the trace of the reconstructed 
solution Note that the iterated quasi-reversibility method gives good result 

even with severely corrupted data. 

Finally, on figure 6, we present the reconstructed Robin coefficient on Fc, which 
was our main objective. Again, the reconstruction is still acceptable for high level 
of noise on the data. 

6.2. One-dimensional heat equation. We now focus on the data-completion 
problem for a one-dimensional heat equation. The problem reads: find u G T) x 
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Figure 3. Residual vs. number of iterations, until stopping crite¬ 
rion is reached, for a = 1%, 2% and 5%. 


dtu 

— ^xx'^ 

in (0, T) X (a, h) 

u{t, a) 

= 9dj 

tG (0,T) 

dxu{t, a) 

= 9n^ 

tG (0,T). 


(a, b)) such that 


Note that, as dxxU = dtu G L^(0, T; F^(a, b)),we havep := dxU € L'^{0,T; H^{a, b)), 
and hence gN{t) = G L^{0,T) without additional assumption, which is not 

the case for the multi-dimensional case. Hence the equivalent data-completion prob¬ 
lem with additional unknown p reads: for igo,g%) G L^(0,T) x L^(0,T), find 
u G T) X (a, 6)) and p G L'^{0,T-, H^{a] b)) such that 


dtu 

= dxP 

in (0,T) X (a, 


dxU 

= P 

in (0, T) X {a. 


u{t, a) 

= 9d^ 

tG(0,T) 


p{t, a) 

= 9n^ 

tG (0,T). 



According to our study, the quasi-reversibility regularization of this problem is: for 
e > 0, find Ug G i?^’^((0,T) x (a, 6)) and pg G L‘^{0,T; H^{a,b)) such that for all 
V G T) X (a, &)), for all g G L‘^{0,T; H^{a, b)) 

j J (^{dtUg - dxPg) {dtv - dxq) + {dxUg - pg) {dxV - q)J dx dt 

+ J (ug{s,a)v{s,a)+Pgia,s)q{a,s)j dt+e J J (jdtUg dtV+dxUg dxV+pg q"j dx dt 
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Figure 4. \us^ — w|, for a = 1%, 2% and 5%. 


= (^9Dis)v{s,a) + g%{s)q{s,a)j ds, 

and the iterated quasi-reversibility method: for e > 0, define = (0,0) 

and for all M € N, x (a, 5)) and € L^(0, T; i/^(a, 6)) are such 

that for all v G iJ^’^((0, T) x (a, 6)), for all q G L‘^{0,T] H^{a, b)), 

j J ({dtu^ - d^p^) {dtv - dxq) -b {dxuf - p^) {d^v - q)j dx dt 

pT r'^ r ^ 

-b J (^u^{s,a)v{s,a)+p^{a,s)q{a,s)^dt+e J j (^dtu^ dtV+dxU^ dxV+p^ q^ dx dt 

= J (^9Dis)v{s,a)+g%(s)q{s,a)^ ds+e J J dtV+dxU^~^ dxV+p^~^ q^ dx dt. 

We discretize the space and L‘^{0,T; H^{a, bj) using a tensorial product of 

Lagrange finite elements, namely 0 P^ finite elements for and x P^ for 

In our simulations, we choose T = 1, a = 1 and b = 2. We consider two 

1 fx^ \ 

exact solution of the heat equation ui(t,x) := - -b a: (1 -b 21) ) and U 2 {t, x) := 


e-*/Sin (t/2). 

The corresponding exact data {go, gN) are corrupted pointwise by a normal noise 
with zero means and variance one, which is scaled so that the noisy data (g^Tdif ) 
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Figure 5. Exact Dirichlet data, noisy Dirichlet data, and trace of 
on r, for a = 1%, 2% and 5%. 


verifies 

\\ 9 d 5-dIIoo Ci|| *?-D IIoo ) \\ 9 n 5Af||oo Q^IISaIIoo- 

In our experiments, we test our method with a = 2% and a = 5%. As in the 
elliptic case, we choose £ = 1, and stop the iterations of the method once the 
stopping criterion is reached. 

In figures 8 and 9, we present the relative error over Q, defined as the ratio 

u^{S) — u 
llltlloo 

for both solutions Ui and U 2 - We see that the iterated quasi-reversibility method 
gives also good reconstruction for this parabolic problem, even for high level of noise 
on both Dirichlet and Neumann data. 

Finally, in figures 10, we present the evolution of the residual quantity 

11^2(2) + -pf 11^2(2) + llwf (■,«) -5Dlli2(0,T) + IIP"(->“) - 5 ^IIl2(0.T) 

during the iterations of the method, until the stopping criterion is reached. 
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